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ABSTRACT 

We present selected topics on heavy-light physics on the lattice in order to illustrate 
the current status of this field. In particular results concerning fg, semileptonic 
decays, the decay B — > K*^f and heavy baryon masses are discussed. Special 
emphasis is paid to the question of systematic effects which stem from the lattice 
discretization of QCD. 

1. Introduction 

The tremendous interest heavy-light quark systems have gained over the recent 
years is mostly due to the need of a precise determination of those Cabibbo Kobayashi 
Maskawa matrix elements, which govern the weak transition of heavy into light 
quarks. In order to extract these fundamental parameters from experimental mea- 
surements, an accurate non-perturbative calculation of the QCD matrix elements 
involved is necessary. The aim of lattice QCD in this context is to compute the re- 
quired QCD part directly from the very definition of the theory itself, i.e. it's action 
Sqcd- 

The non-perturbative "per se" method to deduce QCD properties directly from 
hrst principles would be "simply" to solve the QCD path functional. For example 
the vacuum expectation value of an operator O in terms of this functional is given by 

(0|O(#,*,A)|0) = \- [ D[A]D[V]D[i$]0(y , A)e lS ^ CD[ao ^^ ] . (1) 

Here and g denote the fermion fields, the gluons and the (bare) coupling. As 

it stands, however, Eq. (1) is mathematically ill defined. In order to give it a definite 
meaning one needs regularization. 

The lattice method regularizes Eq. (1) by replacing the space-time continuum 
by an euclidian space-time lattice. In this way everything becomes well defined and 
the path functional can - in principle - be integrated out without any further ap- 
proximation. In actual applications one mostly uses Monte Carlo techniques for this 
purpose. 

*Invited talk given at the 138. WE-Heraeus-Seminar on Heavy Quark Physics, 14. 12-16. 12. 1994, Bad 
Honnef, Germany 
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Of course the lattice method has it's shortcomings. First of all discretization 
means introduction of a finite lattice spacing a, which corresponds to a finite energy 
cutoff a -1 . In order to extract continuum properties from lattice calculations, the 
influence of this cutoff onto the measured quantities must be investigated and finally 
an extrapolation a — > is necessary. Clearly, on a given lattice, masses and energies 
have to stay away from the cutoff. 

Secondly, a finite number of lattice points corresponds to a finite extension of the 
lattice in space and time. This, eventually, induces finite size effects, which also must 
be controlled and removed by an appropriate extrapolation. The great advantage of 
the lattice method is however, that it has in itself the power to control it's systemat- 
ics: Lattice spacing and lattice size can be varied in order to visualize and finally to 
remove the corresponding lattice artefacts. 

The investigation of heavy-light systems with lattice techniques is a very challeng- 
ing task. On the one hand one is forced to keep the physical size of the lattice large 
in order to take care for the light - and thus long ranged - degrees of freedom. On 
the other hand the heavy degrees of freedom require a high lattice resolution, i.e. a 
small lattice spacing. 

Current computer resources enable us to implement lattice resolutions up to a -1 ~ 
2.5 — 3.5 GeV , keeping finite size effects small. Thus, direct calculations in the D 
meson region are feasible today. Beauty properties, however, can be accessed only 
by more elaborate lattice techniques. This " hunting for the beauty" is done with 
different, but interrelated methods. 

The most direct attempt is to investigate the heavy quark mass dependence of the 
required quantities up to approximately two times the charm mass and then to ex- 
trapolate the results to the B meson. This extrapolation can be replaced by a, clearly 
safer, interpolation if one performs in addition the lattice calculation at the heavy 
quark mass limit — > oo, where the heavy quark contribution can be integrated 
out. 

A current, very active line of investigation addresses the question of how to reduce 
lattice artefacts by an improved discretization of the QCD action 1 ' 2 ' 3 . The general 
strategy is to include higher order (in a) terms into the discretized action, which 
weaken the cutoff dependence of the lattice data. Simulations using such an im- 
proved action are still "first principle" calculations, as the additional terms vanish in 
the continuum limit. 

In a less fundamental approach one tries to reduce the lattice artefacts by a suitable 
redefinition of the quark field normalization 4 ' 5 ' 8 . Taking the free case as a guide, one 
successively employs mean field ideas and non-relativistic strategies in order to cor- 
rect for the finite a effects. The benefit of such an ansatz is still under debate, but 
clearly it always can be checked directly on the lattice by variation of a. 

This talk is organized as follows. The next section will report on the status of 
fs on the lattice. Section 3 deals with the semileptonic decays D — > K,K*(vl) and 
B — > 7r, p(vl). The decay B — > K*^ will be discussed in section 4. Section 5 is devoted 
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to heavy baryon spectroscopy. Summary and conclusions will be given in section 6 J 
All the lattice results presented here are obtained in the quenched approximation, 
which - loosely speaking - neglects internal fermion loops. It is used because of 
limited computer resources and - hopefully - it will be removed within the next years. 
Although it's impact cannot be quantified exactly at the moment, it has been proven 
to work accurately at least in the light fermion sector: quenched lattice calculations 6 
of light fermion masses and decay constants are in agreement with the experimental 
results within an error margin of 10%. 

2. Status of /b 

The decay constant fps of a pseudoscalar particle is defined by the QCD matrix 
element 

(0\A o \PS) = fpsm PS , (2) 

with the 0-th component of the (heavy-light) axial current A = ^ /i7o75 v f /• It's 
discretized version (0\A o \B) lat \ which can be evaluated on a given lattice, is related 
to the continuum matrix element by 

(0\A o \PS) = Z A (g 0} m h a } mi a)(U\A \PS) latt . (3) 

The bare coupling g and the quark masses and m\ are the only free parameters 
on the lattice. The lattice spacing a is determined in physical units by fixing one 
lattice quantity to it's experimental value. One frequently uses m p , /„. or <r(string 
tension) to set the scale. 

The renormalization constant Za is composed of a short distance part, which depends 
on the coupling g and can be calculated in perturbation theory, and a long distance 
part, whose impact has to be removed by proper variation of a. 

The question about the size of the systematic error introduced by the perturbative 
calculation of Za is difficult to answer precisely. It has been shown over the recent 
years 7 ' 8 ' 10 that the replacement of the bare coupling g by a suitably chosen effective 
coupling g e jf leads to a clear improvement in the comparison of perturbative results 
and lattice simulations. Non-perturbative studies of Za - including ma effects - are 
also in progress 11 . 

In order to extract the matrix element (0\A o \P S) latt one calculates the 2-point 
correlation function (0|Aj(x, t)A o (0, 0)|0), which can be decomposed into a sum over 
energy eigenstates n with mass m n 

£(0|Aj(f,t)Ao(0,0)|0) =]T|(0|Ao|n)| 2 e -"^. (4) 

The required ground state is obtained by analyzing the correlator at "suitably large" 
euclidian time t 

\{0\A o \PS)\ 2 e - m '" t = lim V \{0\A o \n)\ 2 e- m ^ . (5) 



^We do not discuss here the results of non-relativistic lattice QCD. For a review the e.g. J. Sloan 9 . 
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As the onset of the region of ground state dominance can be ascertained directly by 
analysis of the data, there is no basic problem to extract (0\A o \P S) latt . 

On a finite lattice and with limited statistical accuracy one meets however serious 
difficulties in isolating the ground state unambiguously, especially in the case of heavy 
quarks 12 ' 25 . The cure to this problem comes with the introduction of quark wave 
functions. Use of the exact ground state wave function would yield 

(0|A o |n) = for n ^ groundstate , (6) 

and (0|A o |PS') could be determined already at small values of t. An educated guess 
of the wave function equally leads to significant improvement, reminiscent to the 
variational approach in quantum mechanics. 

2.1. The Static Limit 

Considerable progress has been achieved in the determination of fs in the static 
limit with lattice methods during the last two years. The simulations have been 
improved with respect to their statistical significance - O(f00) gauge configurations 
is standard already - and by the use of quark wave functions, which allow for a much 
more reliable ground state isolation. A variety of lattice results over a reasonable 
range of lattice spacings is available now such that the extrapolation to zero lattice 
spacing is feasible. 

In Fig. 1 we show the collected results of the various groups! First of all we 
observe that the scaling combination Z'^ at f l s ^ t \/Mps exhibits a clear dependence on 
a and an extrapolation to a = is necessary. The PSFWuppertal collaboration 14 
has performed this extrapolation on it's data, with the result f stat = 230 ± 22 ± 26 
MeV. This is nicely consistent with C.R. Allton's 22 result, f stat = 230 ± 35 MeV, who 
included all available data in the analysis. 

A second look at Fig. 1 however reveals, that the data of the FNAL 15 and 
Kentucky 17 groups is significantly lower than the bulk of results at comparable a. 
This observation finds it's explanation in the fact, that these groups have developed 
and applied a new method of projecting onto the low energy eigenstates of Eq. 4. This 
"multistate smearing" method renders a conclusive test of ground state isolation and 
enables for the extraction of the mass gap between ground state and first excited 
state. The a — > extrapolation of the FNAL data yields the result 

fstat = 188 ± 23 ± 15 ± 26 ± 14MeV , (7) 

where the first error is due to statistics, the second gives the scale uncertainty, the 
third is due to the extrapolation in a and the last one estimates the uncertainty 
of Z^ at . Although this result is still consistent with the former analysis mentioned 
above, it's somewhat lower value appears to be the more likely one? 

♦The latest results of C.R. Allton et al. 19 are not discussed here , since a final analysis, including 
systematic uncertainties, is missing yet. 

SR. Sommer showed 23 that the PSI-WUP data indeed moves down to the FNAL results if the 
information about the mass gap is included into the analysis. His re-analysis of the FNAL data, 
however, reveals a much larger uncertainty due to the extrapolation in a. 
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Figure 1: Z% at 'JltltVMps in units of the string te nsion E = era as a function o/E. 
The continuum is at E = 0. The data is taken from: ELC 13 , PSI-WUP 14 , FNAL 15 , 
BLS 4 , UKQCD 16 , Kentucky 17 and APE 18 . In order to compare, all data have been 
uniformly reseated with the string tension taken from Bali and Schilling 20 and with 
the mean field improved 1-loop 21 Z^ at ■ 



2.2. fs with Propagating Quarks 

Finite mass heavy-light systems on the lattice potentially bear the danger of being 
contaminated with large discretization errors, as one is tempted to work near the 
lattice cutoff, where the condition aMps <C f is not valid. 

In order to avoid this contamination one performs the lattice calculation away 
from the cutoff and then extrapolates^the results to the heavy mass region. 

On top of this one tries to suppress the discretization errors by use of improved 
actions or redefinition of the quark held normalization, as explained in the introduc- 
tion. A popular choice 16 ' 18 of such an improved action is the Sheikholeslami-Wohlert 
action 1 , which reduces the discretization effects from O(a) to 0(a/lna) } compared to 
the standard Wilson action 24 . Unfortunately one has to pay for it's application by 
larger statistical noise and an increased need of computer time and memory. Concern- 
ing the redefinition of quark held normalization 4 ' 5 ' 8 the question arises immediately , 
how it can be checked, that such ideas really do improve the situation. Fortunately 
the answer to that can also be given immediately: As improvement means to be - 
at given finite a - closer to the continuum than in the unimproved case, one has to 
study the a dependence of the observables in question. Weaker dependence on a then 
signals improvement. 

As an example of the feasibility of such test we show in Fig. 2 the results of 

^Or interpolates, if the result in the heavy mass limit is available. 
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Figure 2: TTie pseudo-scalar decay constant fps as a function of a, for various me- 
son masses Mps- The circles (squares) refer to the standard (KroMac) normalization. 
The solid lines indicate the extrapolation to the continuum. The continuum result is 
depicted at a a = . 



the PSI-Wuppertal collaboration 25 . They have studied the decay constant fps as a 
function of a and Mps } both in the standard Wilson normalization of quark fields 
and in the KroMac 5 normalization. Obviously the KroMac normalization does not 
weaken the a dependence! 

Finally we present in table f the latest lattice results'lconcerning fp, and fp- It 
is encouraging to see that all results agree within errors, although different methods 
have been used in order to reduce finite a effects. We emphasize that, for the first 
time, it has been possible really to investigate the finite a and - which could not be 
discussed here - the finite size effects. A comparison with the results in the static 
limit (Eq. 7) furthermore reveals, that the large gap between f stat and fp } which had 
been obtained in previous simulations 27 , has almost completely disappeared. This is 
mostly due to (a) the extrapolation in a, (b) the use of proper quark wave functions 
and (c) the increased statistical accuracy of data. 

3. Semileptonic Decays of Pseudoscalar Mesons 

The next step is to calculate the QCD matrix elements relevant for the semilep- 
II For comparison with results from QCD sum rules see e.g. C.A. Dominguez 26 
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Table 1: Latest results for fp, and fs from various groups. a a = oo indicates that an 
extrapolation to the continuum limit has been performed. All groups have calculated 
both f s tat and fps with finite mass quarks. 



group 


fo [MeV] 


fs [MeV] 


method to reduce 
finite a effects 


^[GeV] 


UKQCD 16 




160±g±f| 


SW action 


~ 2.7 


BLS 4 


208(9) ±35 ± 12 


187(10) ±34± 15 


changed 
normalization 


~ 3.1 


PSI-WUP 25 


170(30) 


180(50) 


cont. extrap. 


oo 


MILC 28 


181(4)(18) 


147(6)(23) 


cont. extrap. 


oo 



tonic decays of pseudoscalar mesons. Compared to the calculations of fps this is 
computationally more advanced, as one has to project onto two hadronic ground 
states, which asks for two large time separations on the lattice. Moreover, the need 
of spatial momentum p ^ states leads to an enhanced statistical noise of the lattice 
signals. 

After the early pioneering work of C. Bernard et al. 33 and V. Lubicz et al. 38 , 
much progress concerning the reliability of the results has been achieved over the last 
two years 29 ' 30 ' 31 ' 34 ' 35 . High statistics calculations (O(100) configurations), using quark 
wave functions and fine grained lattices are now available. 

The q 2 dependence of the form factors, which parametrize the required matrix 
elements, has been studied in some detail. For D decays all groups find consistency 
with the pole dominance hypothesis, though within still sizable statistical errors. 

Conventionally, one extrapolates the form factors to their values at q 2 = 0. We 
emphasize that the validity of the pole dominance hypothesis in most cases is not 
crucial for the extrapolation, as the data contains already points near q 2 = 0. The 
corresponding results are collected in table 2** We obtain, that all estimates are in 
reasonable agreement with experiment. 

One might ask whether the existing lattice data allow already for an extrapolation 
a — > 0. In Fig. 3 we investigate this question for the case of f^(0). It is encouraging 
to see that the data, covering already a sizable range in a, exhibit only a weak (if any) 
dependence on the lattice spacing. Due to the large errors however, an extrapolation 
in a appears to be premature. 

In order to obtain the form factors for the decays B — > tt and B — > p } one has to 
extrapolate in the heavy meson mass. In view of the accuracy of the existing data, this 
is a difficult task, but ELC 30 , APE 31 and UKQCD 32 have tackled it. Their findings 
are displayed in table 3, together with a comparison to QCD sum rules. Clearly, the 

*This is a abbreviated version of the nice compilation given by the UKQCD 29 group. A comparison 
with QCD sum rules and potential model calculations is also given there. 
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Table 2: Form Factors at q 2 = for the semileptonic decays D — > K } K*. 



group 






1/(0) 


Ai(0) 


A irt\ 
A 2 (0) 




0.6/_ 8 




1 m+ 30 
l.Ul_ 13 


n 7n+ 7 
0. /0_ 10 


0.66_ 15 


ELC 30 


0.60(15)(7) 




0.86(24) 


0.64(16) 


0.40(28)(4) 


APE 31 


0.72(9) 




1.0(2) 


0.64(11) 


0.46(34) 


BKS 33 


0.90(8)(21) 


0.70(8)(24) 


1 /IQ+45+48 
-L-^O-45-49 


0.83(14)(28) 


n 5Q+14+24 


BG 34 


0.73(5) 


0.73(4) 


1.24(8) 


0.66(3) 


0.42(17) 


WUP 35 


0.76(15) 


0.75(6) 


1.05(33) 


0.59(8) 


0.56(40) 


LMMS 38 


0.63(8) 




0.86(10) 


0.53(3) 


0.19(21) 


Experiment (a) 36 
Experiment (b) 37 


0.77(4) 
0.70(3) 




1.16(16) 


0.61(5) 


0.45(9) 



Table 3: Form Factors at q 2 = for the semileptonic decays B — > 7r, p. The subscripts 
"naive", "HQET" refer to the way in which the extrapolation in raps was done. 



group 


/ + (0) 


V(0) 


MO) 


A 2 (0) 


APE#Q£x 31 
APF • 31 

ELC/jqst 30 
FT f 1 30 

UKQCD 32 


0.29(6) 
0.35(8) 
0.26(12)(4) 
0.30(14)(5) 
24+ - 04 


0.45(22) 
0.53(31) 
0.34(10) 
0.37(11) 


0.29(16) 
0.24(12) 
0.25(6) 
0.22(5) 


0.24(56) 
0.27(80) 
0.38(18)(4) 
0.49(21)(5) 


Ball 39 


0.26(2) 


0.6(2) 


0.5(1) 


0.4(2) 



quality of data has to be increased in order to make accurate predictions. 

4. The Decay B -> A * 7 

To leading order perturbation theory in the weak coupling, the matrix element 
for the decay B — > K*^ is given by 40 

M = e ^0C\(m b )V tb V; s (KVAB) , (8) 

where = scr^q^^ + 75)6, and r\ and q are the polarization and momentum of 
the emitted photon. The non perturbative QCD matrix element (K*\J fl \B) can be 
parametrized by three form factors 41 T;(g 2 ), i = 1,2,3, which need to be evaluated 
at q 2 = 0, since the emitted photon is on-shell. In this limit the form factors obey 
the relations 

T 2 (q 2 = 0) = -iTriq 2 = 0) T 3 {q 2 = 0) = . (9) 
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Figure 3: /^(O) as a function of a. The experimental results are shown at a = 0. 
References are found in Tab. 2. For clearness the data at a a ~ 0.53 GeV~ x have been 
shifted slightly . 

Clearly, the task of lattice QCD is to determine Ti(q 2 = 0) and - for consistency - 
T 2 (q 2 = 0). With current lattice capabilities this requires extrapolation to Mb as well 
as to q 2 = 0. 

As the functional dependence of Ti(q 2 ,Mps) and T 2 (q 2 ,Mps) is not a priori 
known, one has work with ansatze, motivated by pole dominance models and heavy 
quark symmetry. This is not a problem in principle, because the reliability of these 
assumptions can be tested directly on the lattice. Of course the final evidence may 
then depend crucially on the accuracy of data. Three groups 42 ' 43 ' 44 have recently 
studied the decay B — > K*^. The most detailed analysis has been performed by the 
UKQCD collaboration 44 . They find 



0.1591^ ± 0.067 (a) 



T 1 (q 2 = V) 



0.124t^ ± 0.022 (b) 



(10) 



using the assumption of (a) single pole dominance and (b) double pole dominance 
for the q 2 dependence of this form factor. The corresponding value of T 2 (q 2 = 0) is 
consistent with these results if single pole dominance is assumed for T 2 . Unfortunately 
the data is not accurate enough to proof the assumptions. 

The results of Bernard et al. 42 and Abada et al. 43 are consistent (within large 
errors) with the values quoted above. However, the latter ascertain a crucial influence 
of the assumptions made for the q 2 dependence. Thus, much higher accuracy of data is 
needed in order to extract the functional dependence directly from lattice simulations. 

5. Heavy Baryon Spectroscopy 

The determination of baryon masses with lattice methods is a much easier and 
clearer task than the evaluation of form factors. They are extracted from the time 
decay of 2-point correlation functions, which exhibit a comfortably high signal to 
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noise ratio, at spatial momentum p = 0. Moreover one does not have to struggle with 
perturbative renormalization constants, because baryon masses are renormalization 
group invariants. Correspondingly, reliable extrapolations to the b quark mass and 

^ 0.8 
> 

CJ 

''J 0.7 
< 

^ 0.6 
0.5 
0.4 



0.2 0.4 0.6 0.8 1 

1/Mp (GeV -1 ) 

Figure 4: Ma — Mps as a function ofl/Mps and a. The solid line corresponds to a 
global fit assuming no a effects. It gives the value shown with the inverted triangle at 
the B meson mass. The value obtained after extrapolation to the continuum is indi- 
cated by the open circle. f3 = 5.74,6.00,6.26 corresponds to a" 1 ~1 .12, 1.88, 2.78GeV. 

to a = are feasible with today's capacities. In Fig. 4 we display the mass difference 
Ma — Mps as a function of 1/Mps at several lattice spacings a CT , as found by the 
PSI-Wuppertal 45 collaboration. We observe that it behaves nicely linear 1/Mps- The 
finite a effects turn out to be small, which allows for a smooth extrapolation in a. 
They finally get 

M Ab = 5.728 ±0.144 ±0.018 GeV , (11) 
which is consistent with the result of the UKQCD collaboration 46 ^ 

M Ab = 5.900^™ GeV , (12) 

obtained on a single lattice with a a ~ 2.7GeV. 

6. Summary and Conclusions 

In this talk we described the systematic improvements achieved recently on com- 
puting heavy-light physics on the lattice, and presented the latest results from such 
simulations. We have seen that - although the large mass of the beauty quark still 
impedes a direct evaluation - high precision calculations in the range up to two times 
the D mass allow already in some cases for a safe extrapolation to the B meson. 
Due to the advent of fast parallel computers the cutoff dependence of lattice results 

ft A comprehensive list of UKQCD results concerning beautiful Baryon masses is given also in that 
reference. 




a = 6.26 
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can be studied now in detail, which makes the extrapolation to the continuum fea- 
sible. Given this necessary condition, the door is open to all kinds of improvement 
ideas, whose benefit can be tested with the "lattice lab". 

Current lattice capabilities have enabled us to consolidate the value fs = 180(50)MeV 
and to give reliable estimates for the form factors of semileptonic D decays. First 
high statistics results concerning the decays B — > 7r, p } B — > K*^ and the masses of 
beautiful baryons are also available now. 

One might wonder, whether the "hunting for the beauty" will end up in a "ren- 
dezvous". In the worst scenario, the hunter will have to wait three years, to the 
advent of 60 Gflops machines, which are powerful enough to increase lattice resolu- 
tion from a -1 ~ 3GeV to a -1 ~ 6GeV. But there are many good ideas to reach the 
goal in a shorter time. 
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